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Abstract 

In models with low-energy supersymmetry breaking, it is well-known that charged slep- 
tons can be significantly lighter than the lightest neutralino, with the gravitino and lighter 
stau being the lightest and next-to-lightest supersymmetric particles respectively. We give 
analytical formulas for the three-body decays of right-handed selectrons and smuons into 
final states involving a tau, a stau, and an electron or muon, which are relevant in this 
scenario. We find that the three-body decays dominate over much of the parameter space, 
but the two-body decays into a lepton and a gravitino can compete if the three-body phase 
space is small and the supersymmctry-brcaking scale (governing the two-body channel) is 
fairly low. We study this situation quantitatively for typical gauge-mediated supersymme- 
try breaking model parameters. The three-body decay lengths are possibly macroscopic, 
leading to new unusual signals. We also analyze the final-state energy distributions, and 
briefly assess the prospects for detecting these decays at CERN LEP2 and other colliders. 
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1. Introduction 



Supersymmetry-breaking effects in the Minimal Supersymmetric Standard Model (MSSM) 
are usually introduced explicitly as soft terms in the lagrangian. In a more complete theory, 
supersymmetry is expected to be an exact local symmetry of the lagrangian which is sponta- 
neously broken in the vacuum state in a sector of particles distinct from the MSSM. There are 
two main proposals for how supersymmetry breaking is communicated to the MSSM particles. 
Historically, the more popular approach has been that supersymmetry breaking occurs at a scale 
> 10 10 GeV and is communicated to the MSSM dominantly by gravitational interactions. In 
this case, the lightest supersymmetric particle (LSP) is naturally the lightest neutralino (Ni). 
One of the virtues of this gravity-mediated supersymmetry breaking scenario is that a neu- 
tralino LSP can easily have the correct relic density to make up the cold dark matter with a 
cosmologically acceptable density. 

Recently, there has been a resurgence of interest in the idea [H [2] that supersymmetry- 
breaking effects are communicated to the MSSM by the ordinary SU(3)c x SU(2)l x U(1)y 
gauge interactions rather than gravity. This gauge-mediated supersymmetry breaking (GMSB) 
scenario allows the ultimate supersymmetry-breaking order parameter y/~F to be much smaller 
than 10 10 GeV, perhaps even as small as 10 4 GeV or so, with the important implication that the 
gravitino (G) is the LSP. The spin-3/2 gravitino absorbs the would-be goldstino of supersym- 
metry breaking as its longitudinal (helicity ±1/2) components by the super-Higgs mechanism, 
obtaining a mass = F/^/3M P = 2.37(^/100 TeV) 2 eV, where M P = (8vrG N ewton)" 1/2 = 
2.4 x 10 18 GeV is the reduced Planck mass. The gravitino inherits the non-gravitational interac- 
tions of the goldstino it has absorbed [3j- This means that the next-to-lightest supersymmetric 
particle (NLSP) can decay into its standard model partner and a gravitino with a characteristic 
decay length which can be less than of order 100 microns (for y/F < 10 5 GeV) or more than 
a kilometer (for y/F > 10 7 GeV), or anything in between. This leads to many intriguing phe- 
nomenological possibilities which are unique to models of low-energy supersymmetry breaking 
[3-10]. For kinematical purposes, the gravitino is essentially massless. The perhaps surprising 
relevance of a light gravitino for collider physics can be traced to the fact that the interactions 
of the longitudinal components of the gravitino are the same as those of the goldstino it has 
absorbed, and are proportional to I/tHq (or equivalently to 1/F) in the light gravitino (small 
F) limit [3]. 

In a large class of models with low-energy supersymmetry breaking, the NLSP will either 
be the lightest neutralino or the lightest stau (fi) mass eigenstate. Our convention for the stau 
mixing angle 9f is such that 

( h \ = ( cos 6 f sin 0f\ ff L \ , . 

\f 2 ) ~ \-sm9 f cosOf) \f R ) 1 > 
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with mfj < mf 2 and < Of < ir (so sin Of > 0). The sign of cos T depends on the sign of /i (the 
superpotential Higgs mass parameter) through the off-diagonal term —[im T tan f3 in the stau 
(mass) 2 matrix. This term typically dominates over the contribution from the soft trilinear 
scalar couplings in GMSB models, because the latter are very small at the messenger scale 
and because the effects of renormalization group running are usually not very large. For this 
reason, it is quite unlikely that cancellation can lead to cos Of ~ in these models, unless the 
scale of supersymmetry breaking is quite high. In GMSB models like those in Ref. [9] which 
are relevant to the decays studied in this paper, | cos Of | ranges from about 0.1 to 0.3 when the 
mass splittings between f\ and the lighter selectron and smuon are less than about 10 GeV. 
That is the situation we will be interested in here. The selectrons and smuons also mix exactly 
analogously to Eq. ([1]). However, at least in GMSB models, their mixings are generally much 
smaller, with cos Op,/ cos^f ~ U^/Ut ~ 0.06 and cos^g/ cos6*f ~ Ve/Vr ~ 3 x 10 -4 . Therefore, 
in most cases one can just treat the lighter selectron and smuon mass eigenstates as nearly 
unmixed and degenerate states. We will write these mass eigenstates as eR and fiR, despite 
their small mixing. We will also assume, as is the case in minimal GMSB models, that there 
are no lepton flavor violating couplings or mixings. 

The termination of superpartner decay chains depends crucially on the differences between 
m^, mfj, and rug (in this paper £ is generic notation for e or /i). We assume that -R-parity 
violation is absent, so that there are no competing decays for the NLSP. If the NLSP is N\ with 
JTijy- < Tifi ~ m r, then the decay N\ — > 7G can lead to new discovery signals for supersymmetry, 
as explored in Refs. [3-9]. In other models, one finds that the NLSP is t\ [6]. Here one must 
distinguish between several qualitatively distinct scenarios. If tan j3 is not too large, then eR 
and £lr will not be much heavier than fi, and the decays £r — > lrf\ and £r — ► £N\ will not 
be kinematically open. In this "slepton co-NLSP scenario", each of en, (Ir, and fi may decay 
according to e.R — ► eG, jiR — > fiG and fi — > tG, possibly with very long lifetimes. There can 
also be competing three-body decays Ir — > V(p r f\ through off-shell charginos (Ci). However, 
these decays are strongly suppressed by phase space and because the coupling of Ir to v^Ci is 
very small. In the approximations that m? — m? <C m% and 1 — m? /m~ <C 1, one finds 

V Tl '1 1,1 '1 V Tl 



f = £,r. Here Uij is one of the chargino mixing matrices in the notation of [11] and g is 
the SU(2)l gauge coupling. (Of course, the decay i\ — > ugv T fi has the same width.) For JIr 
decays, we find that this width is always less than about 10 -7 eV in GMSB models like the ones 
discussed in [9] if mp a — mf x < m T and mp R J> 80 GeV. The maximum width decreases with 
increasing mp R as long as we continue to require that the decay \1r — * fJ,rf% is not kinematically 
open. (For the corresponding cr decays, the width is more than four orders of magnitude 
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UncosOf— (yf/g)Ui2 sin 0?, with Yukawa couplings yj = gmj /(^/2my/ cos/3), for 
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Figure 1: Right-handed selectrons and smuons can decay according to £ R — > £ T + f x or £ R —* 
£~T~f^, with different matrix elements, through virtual neutralinos iVj (i = 1,2,3,4). 

smaller.) This corresponds to physical decay lengths of (at least) a few meters unless the 
sleptons are produced very close to threshold. It is possible to have somewhat enhanced widths 
if tUq — is decreased or if cos 9^ is increased compared to the values typically found in 
GMSB models. However, even if the decays £r — > V(p T T\ can occur within a detector, they will 
be extraordinarily hard to detect because the neutrinos are unobserved and the f\ momentum 
in the lab frame will not be very different from that of the decaying £r. The subsequent decays 
f\ — * tG can be distinguished from the direct £r — ► £G, but if the latter can occur within the 
detector, then they will likely dominate over £r — » v£P T fi anyway. So it is very doubtful that 
the decays £r — > V(p T f\ can play a role in collider phenomenology. 

For larger values of tan /3, enhanced stau mixing renders f\ lighter than e# and fin by 
more than m T . In this "stau NLSP scenario" , all supersymmetric decay chains should (naively) 
terminate in f\ — * tG [6j [TOl E] , again possibly with a very long lifetime^ If the mass ordering 
is niji R — and/or mg R — m e > , then the two-body decays (1r — > [iN\ and/or en — ► eN\ 
will be open and will dominate. In the rest of this paper, we will instead consider the situation 
in the stau NLSP scenario in which m^> m| — mi > + m T . In that case, {xr and/or en 
can decay through off-shell neutralinos in three-body modes Jar — > \xrf\ and/or eft erfi, as 
shown in Fig. [TJ Here one must be careful to distinguish between the different charge channels 
T + f^ and T~f^ in the final state, for a given charge of the decaying slepton. In the following 
we will give formulas for T(£ R — > £~t + t±) and T(£ R — ► £~t~t^~), which in general can be quite 
different^ except when the virtual neutralino is nearly on shell. [Of course, these are equal to 
r(f^, — ► £ + t~t± ) and T(£^ — > £ + t + t 1 _ ), respectively.] These three-body slepton decays have 
been rightly ignored in previous phenomenological discussions of the MSSM with a neutralino 
LSP, in which the two-body decays £r — ► IN\ (and possibly others) are always open. However, 
in models with a gravitino LSP, N\ is allowed to be much heavier, so it is important to realize 
that three-body decays of en and fin are relevant and can in principle imply long lifetimes and 

An important exception occurs if \m,f 1 — m^J < m T and < m; . In this "neutralino-stau co-NLSP 

scenario" , both of the decays fi — * tG and Ni —> 7G occur without significant competition. 
2 We are indebted to Nima Arkani-Hamed for pointing this out to us. 
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macroscopic decay lengths. In the following, we will present analytical results for the three-body 
decay widths of eR and and study numerical results for typical relevant model parameters. 



2. Three-body slepton decay widths 



Let us first consider the "slepton-charge preserving" decay £ R — > t t+t-l , keeping t mixing 
effects. The matrix element for the relevant Feynman diagrams in Fig. 1 can be written as 

• M = E <Pi)(Pn4 + PLbT) lm2 J* Jl r2)1 (M + PR^vipr), (3) 
where P L)R = (1 ± 7s)A and x t = 2pz ■ p e /pj , and 



"J 



J = V2g' N n sin 6? + y T N j3 cos 6 f ; (4) 
1 

V2" 



b] = - (g'N^ + gN^) cos 6 f + y T N* 3 sm9f, (5) 



with exactly analogous formulas for and bj, with f £. Here we have adopted the notation 
of Ref. [11] for the unitary (complex) neutralino mixing matrix iVjj with all real and 

positive, and g and g' are the SU(2)l and U(l)y gauge couplings. Our fermion propagator is 
proportional to (— p 1 + m)/(p 2 + m 2 ), with a spacetime metric signature (—+++). 

Summing over final state spins and performing the phase space integration, we obtainj§ 

rt^^l^EE^ (6) 

i,j=l a=l 



in terms of coefficients 



cf)=J*a*44* + bfb*libl\ (7) 

$=4b]ajbr + bfapiar, (8) 

eg) = 2Re[af fejafaf* + &f ojftfftf ], (9) 

cf) = -2Re[af 6J&f&f* + bf ojafof*], (10) 

cg ) = -4Re[6faJaf6r], (H) 
cg } = -4Re[af aj&f&f*], (12) 



3 Similar formulas can be derived for the three-body decay widths of all sfermions in the MSSM. Here we have 
neglected higher order effects, including contributions to the neutralino widths from final states other than tti, 
since we will be interested in the situation in which is not too close to Ttif . 
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and dimensionless integrals I-?' defined as follows. First, we introduce the mass ratios rf = 
m^/mg^, r T = m T /m^ R , re = me/rug , and r^= m ft. /rug with ri <C r T <C rf < 1 — r T — re 
and rft > l — r£. Then we find 

iff = J dx e (x e - 2r|)(l - x t + rf )(1 -x t + rf + r 2 T - r 2 f )f ih (13) 

4 2) = r iv/iv, / dx t " 2rf)(l - x e + rj + r 2 T - r?)/y, (14) 



4 3) = 2r T r^. y da* (x< - 2r?)(l - x e + rj)f t] , (15) 

4 4) = 2r t r* f dx e (1 - x e + rf )(1 -x e + rj + r 2 T - r^fa, (16) 



l[f = 2r e r T r^r^ I dx e (l-xe + rj)fij, (17) 

(6) 



2r^r T y da* (1 - a* + rf )%•, (18) 



where 



^-4r|AV2[i-^ + 



r 2 ,r 2 ,r|] 



(1 - x e + rj) 2 (rjL_ - 1 + X£ - r|) (r|. - 1 + x £ - rj) 



(19) 



with A 1//2 [a, 6, c] = Va 2 + fr 2 + c 2 — 2ab — 2ac — 2bc. The limits of integration for xe are 2re < 
xe < 1 + rj — r 2 — r 2 — 2r T rf. The matrix element and decay width for the "slepton-charge 
nipping" channel £7, — > l~r~ft are obtained by replacing — ► 6^* and 6^ — ► a^"* everywhere 
in the above equations. 

In GMSB models like those studied in Ref. [9] which are relevant to these decays, one finds 
that rrig R — m^ R is at the most a few tens of MeV, so we will neglect the distinction between 
rrig R and mn R . It is an excellent approximation to take = except when the mass difference 

Am = mji — rrify — m T (20) 

is a few hundred MeV or less, and r e = is of course nearly always a good approximation. It is 
also generally an excellent approximation to neglect smuon and selectron mixing and Yukawa 
couplings in the matrix element, so that V2g'Nj± and b] » 0@ The effects of ij 4,5 ' 6 ^ 

are usually quite negligible because of the re and tf- suppressions. An instructive limit which 
is often approximately realized in GMSB models (or, in generic models with gaugino mass 
unification, whenever is sizeably larger than the gaugino mass parameters) is the case in 
which the contributions from a Bino-like N\ dominate, with ~ 0.5m^ 2 <C m^, m^. Since 
the decaying £r essentially couples only to the Bino (B) component of the virtual neutralinos, 
this approximation is quite good for a large class of models where \N\x\ is not too far from 1. In 



4 We have calculated the effect of including the smuon mixing and the muon Yukawa to be at the level of a 
few to ten percent of the total smuon width, for typical GMSB models from Ref. 
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that case, we may neglect the contributions of virtual N% and N± because of the coupling 
constant suppressions together with the suppressions due to larger neutralino masses. With 
these approximations, the expressions for the decay widths simplify to 



and 



R 



-T+ff) 



512tH 



AA 2 !^ + IS1P/. 



2,(2) 
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2Re(A 1 Bt)l[f 



A 1 = 2g' 2 \N n \ 2 sinflf + V2g'y T N^N 13 cos6 f , 

Bi = g' 2 N{l cos Of + gg'N^Nfz cos 9 f - V2g'y T N^N^ sin Of, 



-T-f+) 



B[\ 2 I$ + K| 2 /{^ - 2Re(A' 1 B' 1 *)l£ 



|2 r (2) 



• ur- 5127r s 
A; = 2g' 2 Nf 1 sin 6» f + V2g'y T N u N 13 cos f , 

= ^iJVn | 2 cos Of + gg'N n Nl 2 cos 0? - V2g'y T N n N? 3 sme?. 



(3) 



(21) 

(22) 
(23) 

(24) 

(25) 
(26) 



We will be interested in the situation in which Am is small (less than 10 GeV). This implies 
that tan/3 is not too large0 and thus f\ has a large tr content. However, as we will see in 
the next section, it is usually not a good approximation to neglect stau mixing altogether (by 
setting sin^f = 1, cos#f = everywhere), because | cos#f | is likely to be at least 0.1 as we have 
already mentioned. Near threshold, the range of integration includes only small values of xi, 
so that the dimensionless integrals /q and /j^ and scale approximately proportional to 
l/( r ^ ~~ I) 2 an d r \ /( r % ~ I) 2 an d r N 1 /( r % ~ I) 2 ' respectively This means that the decay 
width is suppressed as (or equivalently m^) is increased, with other parameters held fixed. 
This is particularly likely in GMSB models with a large messenger sector and a high scale of 

(2) (3) 

supersymmetry breaking. Furthermore, the relative sizes of the / xl and / xl contributions are 
enhanced in the large limit. It is important to note that as is increased, T(£^ — ► £~T~f^ ) 
becomes larger than T(£^ — ► £~r + fi), because of this effect together with the fact that A± 
and A'i typically have larger magnitudes than B\ and B 1 . Note also that the / xl contributions 
appear to be suppressed by a factor of r T , but this turns out to be illusory since near threshold 
m T is not the only small mass scale in the problem; in particular it can be comparable to or 
even much larger than Am — m^ which determines the kinematic suppression of the decay. 



3. Numerical results 



Some typical results are shown in Figs. [2][5l In Fig. [21 we give the total three-body decay 



5 For example in the GMSB models studied in [9j with < Am < 10 GeV, the relevant range for tan /3 is from 
about 5 to 20 for sleptons that could be accessible at LEP2 or Tevatron upgrades, with smaller values of tan j3 
corresponding to smaller Am. 



7 



10 4 
10 3 
10 2 

1 10" 

<3io -1 
u 

10' 2 
10" 3 
10" 4 



Figure 2: The decay widths in meV for e.R — * erf\ (solid lines) and /1r — > fj,rf\ (dashed lines), 
including both T + f-f and T~fi final states, as functions of Am = — — m T . The widths 
have been computed using Eqs. (f2Tj) - (f26j) with = 90 GeV, cos8f = 0.15 and = 1.1, 1.5, 
2.0, and 3.0 (from top to bottom), with the approximation Nn = 1 and N12 = N13 = 0. 

widths for e.R and fin (including both t + t^ and t~T\ final states) as functions of Am for 
= 90 GeV and four choices = 1.1, 1.5, 2.0, and 3.0. (In the GMSB models studied 
in Ref. [9], one finds ^ 1.8, but it is possible to imagine more general models with 
larger values.) Here we have chosen the approximation of Eqs. (|2ip -(|26p with Nn = 1 and 
cos6» f = 0.15. We use m T = 1.777 GeV, = 0.1057 GeV, a = 1/128.0, and sin 2 9 W = 0.2315. 
Realistic model parameters can introduce a significant variation in the decay widths, and in 
general one should use the full formulas given above for any specific model. Our choice of a 
positive value for cos6*f in this example leads to a suppression in the width compared to the 
opposite choice, because of the sign of the interference terms proportional to 7 X1 in Eqs. (|2ip 
and (|24p . These interference terms are often of the order of tens of percent of the total width, 
showing the importance of keeping the stau mixing effects if real accuracy is needed. 

The important ratio of the partial widths for the two charge channels r(^ — » £~t~t^ ) /T(^ 
t~T + f^) is shown in Fig. [3]for the case I = e, as a function of . Here we have chosen values 
of cos#f = —0.3, —0.1, 0.1 and 0.3, and other parameters as in Fig. [2J As expected, this ratio 
is close to 1 when the virtual neutralino is nearly on-shell, and increases with . It scales 
roughly like 7^ /-^ji ~ r \ , U P to significant corrections from the interference term(s). This 
increase tends to be more pronounced for larger cos#f in these models. Because large also 
corresponds to longer lifetimes, the decay — > l~T~f^ is likely to dominate if the three-body 
decay lengths are macroscopic. 
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Figure 3: The ratio T(e^ — ► e~r~f^ )/r(e^ — > e'r+fj") is shown as a function of , for 
four values of cos#f = 0.3, 0.1, —0.1 and —0.3 (from top to bottom). The widths have been 
computed using Eqs. ([2T|) -([26 |) with mg fl = 90 GeV, and with the approximation Nu = 1 and 

#12 = #13 = 0. 

The variation with the stau mixing angle is further illustrated in Fig. HJ where we show the 
total three-body decay width T{e.R — > erfi) including both charge final states with rrie R = 90 
GeV and Am = 1.0 GeV, r$ = 1.1, 1.5, 2.0, 3.0, for the range -0.5 < cos# f < 0.5. Note that 
the total width can vary by a factor of two or more over this range. Here it should be kept in 
mind that at least in the GMSB models studied in Ref. [9], one finds 0.1 < | cos#f | < 0.3, so 
that the whole range shown may not be relevant. In Fig. we show contours of constant total 
three-body decay widths T(e^ — ► erfi) in the Am vs. mg R plane, for the choice = 1.5 and 
cos6*f = 0.15. In both figures we continue to use N\\ = 1, N\2 = #3 = in Eqs. (f2"Tj) - (|2f)j) . 
However, it should be emphasized that in realistic models the effects of deviations from this 
simplistic approximation can be quite appreciable, especially since |-/V"n| 2 can easily be of order 
0.7 or somewhat less in GMSB models, and the width scales essentially like |iVii| 4 . 

As can be seen from these figures, the physical three-body decay lengths for en and JIr 
can be quite large if Am is less than a few GeV and/or m^/m^ is large. In the lab frame, 
the probability that a slepton £r with energy E will travel a distance x before decaying is 
P{x) = e~ x l L ', where 

/ r \ _1 / e 2 \i/2 

i = a2 (r^j mm - (27) 

For sleptons pair-produced at LEP2 (or at a next-generation lepton collider), E in Eq. (|27p is 
simply the beam energy. So if Am is less than a GeV or so (depending on and the specific 
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Figure 4: The dependence of T(eR — ► erfi) (including both charge final states) on cos#f, 
computed as in Fig. [2] but with Am = mg R — mf — m T held fixed at 1.0 GeV. The four curves 
correspond to = 1.1, 1.5, 2.0, and 3.0 (from top to bottom). 



couplings of the model), en and fiR could have a macroscopic and measurable decay length. If 
Am is of order 100 MeV or less, the decay length could even exceed the dimensions of typical 
detectors. 

It is also important to realize that the dominant decay for (.r is not a priori known, since 
the three-body decays £r — > lrf\ have to compete with the two-body decays to the gravitino 
£r — > IG. The latter have a width given by 

T(£ R -> IG) = m\j\§itF 2 . (28) 

For a given set of weak-scale MSSM parameters leading to a calculable three-body width for 
£r, the two-body width Eq. ([28]) is essentially an independent parameter, depending on y/F 
(or on the gravitino mass in "no-scale" supergravity models [13]). For example, for the sets of 
parameters and corresponding widths in Fig. 2, the three-body decay dominates for y/F > 10 3 
TeV for Am — mi down to a few hundred MeV. Alternatively, the minimum possible value of 
yF of order 10 TeV in GMSB models corresponds to a maximum width for £r ^ IG of order 20 
eV (for mi R of order 100 GeV), so Am is expected to be larger than of order 10 GeV before the 
three-body decay dominates. In many of the GMSB models that have actually been constructed 
including the super symmetry-breaking sector [21 [12], this limit is not saturated and yF is orders 
of magnitude larger than 10 TeV, so the three-body decay is expected to dominate unless the 
mass difference is correspondingly smaller. Conversely, in "no-scale" models, the two-body 
decay width might even be much larger than the tens of eV range. 
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Figure 5: Contours of constant total decay width T(cr — > erfi) (from left to right, 0.0001, 
0.001, 0.01, 0.1, 1, 10, 100 and 1000 meV), including both charge channels for the final state, 
and computed with = 1.5 and cos9f = 0.15 with the same approximations as in Figs. [2][D 

4. Energy distributions 



If the three-body decays of e# and Jar indeed dominate, then the £ and r emitted in the 
decay can be quite soft if Am is small. Hence, it is important to address the lepton detectability 
and, in general, the ability to recognize a three-body decay pattern in a real experimental 
environment. Using CompHEP 3.2 |15j plus an implementation of the MSSM lagrangian [16J, 
we have examineqj the (s)particle energy distributions; those of e or /i and r are shown in 
Fig. E](a) and[6tb). Here, we have plotted the results for £~^ — > £~t~t + , but we have checked 
that the shapes of the normalized distributions for 1^ —* £~t + t~ are essentially identical. First, 
we consider a model with N\ = B, = 90 GeV, = 1.1, cos6*f = 0.15, as in the first case 
of Fig. [2j with Am = 1 GeV. Fig. [6ja) shows that the final e or \x (solid thick or dashed line) 
usually has an energy greater than half a GeV in the rest frame of the decaying selectron or 
smuon. Hence, especially when Ir is produced near threshold (as could happen, e.g., at LEP2) 
and the boost to the lab frame is small, a successful search for the e or [i in this model requires 
a detector sensitivity at the level of 1 GeV or better (with low associated energy cuts). The 
r (circles and dot-dashed line) gets most of the remaining available energy, so that E T — m T 
is usually less than 0.5 GeV, while the momentum \p T \ is usually < 1.5 GeV in the £r rest 
frame. It is interesting to note that the final f\ can get up to only 2 GeV in momentum (and 

6 Note that we have checked in great detail and for a wide range of parameters that the partial widths for 
three-body £r decays obtained with CompHEP are in excellent agreement with our analytical results given above. 
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Figure 6: Lepton energy distributions in the rest frame of the decaying to 1~t~t£. (a) 
Normalized distributions for both the final i (solid and dashed lines) and r (circles and dot- 
dashed line) for an ideal model with N\\ = 1; N±2 = N13 = 0; rug = 90 GeV; Am = 
m t R ~ mf i ~~ m " r = ^ GeV; cos#f = 0.15 and = 1.1. Distributions for the other charge 
channel are almost identical. The solid line and the circles (dashed and dot-dashed lines) refer 
to the case £ = e [l = //). (b) The logarithmic version of the solid thick curve in (a) compared 
to normalized electron-energy distributions in four GMSB models chosen from Ref. |9j (thin 
lines). Am is 0.16, 0.30, 2.2, 9.7 GeV respectively from left to right, other details can be found 
in the text. 

usually less), in this case. In the particular model we are considering here, L is of order 5/xm 
at LEP2 [from Eq. (|27p ], and so the kink is impossible to detect. However, the decay length 
could easily be longer in models with, for example, a larger ratio m^/m^ R with fixed external 
particle masses. In those cases where the final leptons are too soft to be detected, the presence 
of such a kink in the charged track might still signal a three-body decay pattern. 

Most of the above considerations strictly apply to the particular model we are considering 
with Am = 1 GeV. Since the prospects for detection depend crucially on Am, it is important 
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to understand how the distributions scale while varying Am (and also other parameters). We 
find that the shapes of the energy distributions in Fig. [6](a) stay basically the same when Am 
is changed, after performing a suitable rescaling of the axes. In addition, we have checked that 
they are only slightly affected by, e.g., changes in and/or stau mixing angle (within models 
with Ni rj B). Only when m^ gets very close to and/or | cos#f | J> 0.3 can deviations 
exceed a few percent (larger deviations are often in the direction of shifting the maximum of 
the e or /u distribution towards slightly lower values, and vice- versa for the tau distribution). 

More generally, in Fig. [U[b), we illustrate the scaling using particular GMSB models from 
Ref. [9] that are relevant for the slepton three-body decays. We show the logarithmic and nor- 
malized electron energy distributions for four models (thin lines) compared to that of Fig. [6|a) 
(thick line). These four GMSB models have, respectively from left to right: = 75.8, 89.8, 
63.7, 69.7 GeV; Am = 0.16, 0.30, 2.2, 9.7 GeV; cos6 f = 0.13; 0.12, -0.21, -0.31; m^ = 97.9, 
95.0, 64.6, 75.1 GeV; |iV ia | 2 = 0.88, 0.97, 0.50, 0.73; so that T(e^ e~T~f^) = 9.11 x lO" 6 , 
1.11 xl0~ 3 , 6.05, and 237 eV and T(e^ -> e~r+ff ) = 5.19 xl0~ 6 , 9.75 xlO" 4 , 5.47 and 171 eV 
[using Eq. © and the corresponding equation for T(f^ — > £~T~f±)]. They were picked in such 
a way as to probe various regions of the GMSB parameter space allowed for models within reach 
of LEP2. Fig. E^b) shows that, in addition to slight deformations of the shapes of distributions 
due to small — 1 and/or large | cosft? | > 0.3, values of |iVii| 2 < 0.7 can produce further 
small changes (as evident from the two models more on the right with larger Am). The total 
deviations are, however, still small enough to allow a model-independent generalization of the 
discussion above concerning the detectability of the three-body decay. Thus, it is expected that 
in most models the e or \i will typically get more than half of the available energy, and hence 
the chance for detection increases straightforwardly with increasing Am. However, the decay 
length of the e# or fiR will drop in correspondence with the total width increase, diminishing 
the chance of detecting a kink in the charged track. Alternatively, for smaller Am, detection 
of the e or /j, (and also the £r kink) is more difficult, but of course the decay length is longer, 
increasing the chance that a kink can be seen. 



5. Discussion 



At LEP2, the process e + e~ — ► 7^ f{~ is the most kinematically-favored one for super- 
symmetry discovery in the stau NLSP scenario. If the decay f\ — ► tG takes place outside 
the detector (or inside the detector but with a decay length longer than a few cm), then the 
stau tracks (or decay kinks) may be directly identified [H [H]. If en and jln can also be 
pair-produced, then the decays £r — > trr\ studied here can come into play, leading to addi- 
tional events e + e~ — > l + 1~ t + t + f7~ or 1 + £~t~t~t^t^ or £ + £~t + t~t 1 i tT. Note that when 
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— > ^ _ t~Ti~) is larger than — > £~T + f{"), the same-sign T ± r ± f{ F f{ F signals are sup- 
pressed compared to the opposite sign signals r + r _ f^"fj _ . In Ref. [9], it was observed that the 
f^R^R production cross section in these models is often significantly larger than that for e^e^, 
because of the interference effects of a heavier neutralino in the t-channel diagrams contributing 
to the latter process. Therefore, one may expect more /z + ' \i~ ~ttt\t\ events than e + e~ ttt\T\ 
events, although this is not guaranteed. We have seen that if Am is smaller than order 1 GeV, 
then the identification of soft leptons and taus may be challenging. However, we noted that in 
just this case the decay length of £r may well be macroscopic, leading to another avenue for 
discovery. Also, since £r decays isotropically in the rest frame, and pair-produced sleptons gen- 
erally do not have a considerable preference for the beam direction, we expect the probability 
for the final particles to be lost down the beam pipe to be small. This is especially true for 
I = fA, where the production does not receive contributions from i-channel neutralino exchange 
(see, e.g., Ref. [9]). 

If fx decays to tG with a decay length shorter than a few cm, then f\ decay kinks will 
be difficult to observe directly at LEP2. Instead, f± r-f production leads only to a signal 
T + T~lfi. This has a large background from W + W~ production, but it may be possible to defeat 
the backgrounds with polar angle cuts [9]. If £r pair production is accessible and £r — ► £G 
dominates over £r — > £tt\, then the model will behave essentially like a slepton co-NLSP 
model, even though the mass ordering is naively that of a stau NLSP model. We have seen 
that this might occur even for a multi-GeV Am. Then the most likely discovery process may 
be e + e~ — > Ah/^r — ► n + n~ as discussed in Ref. [9j. If the decay t\ — > tG is prompt but the 
decays £r — > £rf\ discussed here still manage to dominate over £r — ► £G, then one can have 
events e+e" -> £\£ R t + t + (£ + £-t~t~)$ or t~t~ (£ + Tt + t + )$ or t+t-{£ + £-t + t-)$, with 
the leptons in parentheses being much softer. The first two should have very small backgrounds, 
as will the last one if the soft leptons are seen. 

At the Fermilab Tevatron collider, sleptons can be pair-produced directly or produced in 
the decays of charginos and neutralinos. If the decays f\ — > tG and £r — ► lrf\ both take place 
over macroscopic lengths, then pp — ► C\C\ or C\Ni can lead to events with leptons + jets + 
heavy charged particle tracks (possibly with decay kinks). It is important to realize that both 
the production cross-section and the detection efficiency for such events will likely be greater 
than for the direct production processes pp — > £r£r and f\T\. If f\ — > tG has a macroscopic 
decay length but the decays £r — ► Itt\ studied here are prompt, then there will be some events 
with extra soft leptons and taus. However, the latter may be difficult to detect, and furthermore 
one may expect that C\ and Ni will decay preferentially to t\v t and f\T (or v t t and tv^v) 
rather than through £r. Similar statements apply for the CERN Large Hadron Collider, except 
that the most important source of sleptons may well be from cascade decays of gluinos and 
squarks; in some circumstances those decays may be more likely to contain £r channels. 
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In this paper we have studied the three-body decays of selectrons and smuons in the case 
that the neutralino is heavier. In GMSB models and other models with a gravitino LSP, 
these decays may play a key role in collider phenomenology. In particular, we found that the 
corresponding decay lengths may be macroscopic and the competition with the decays £r — ► IG 
may be non-trivial. We also found that the electron or muon in the final state of the three body 
decay usually carries more than half of the available energy in the rest frame of the decaying 
slepton. 
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